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^S) ' Abstract. The equations of the three-dimensional viscous, compressible, and heat 

conducting magnotohydrodynamic flows are considered in a bounded domain. The 

CS , . viscosity coeflicients and heat conductivity can depend on the temperature. A solu- 

^^ . tion to the initial-boundary value problem is constructed through an approximation 

^^ ' scheme and a weak convergence method. The existence of a global variational weak 

O^ ' solution to the three-dimensional full magnetohydrodynamic equations with large data 

^sj ^ is established. 

< 

2 ' 1- Introduction 

Magnetohydrodynamics. or MHD, studies the dynamics of electrically conducting fluids 
and the theory of the macroscopic interaction of electrically conducting fluids with a mag- 
netic field. The applications of magnetohydrodynamics cover a very wide range of physical 
areas from liquid metals to cosmic plasmas, for example, the intensely heated and ionized 
fiuids in an electromagnetic field in astrophysics, geophysics, high-speed aerodynamics, 
lO ' and plasma physics. Astrophysical problems include solar structure, especially in the outer 

\l I layers, the solar wind bathing the earth and other planets, and interstellar magnetic fields. 

^^ ■ The primary geophysical problem is planetary magnetism, produced by currents deep in 

^D I the planet, a problem that has not been solved to any degree of satisfaction. Magnetohy- 

QP ■ drodynamics is of importance in connection with many engineering problems as well, such 

as sustained plasma confinement for controlled thermonuclear fusion, liquid-metal cooling 
of nuclear reactors, magnetohydrodynamic power generation, electro-magnetic casting of 
K^ \ metals, and plasma accelerators for ion thrusters for spacecraft propulsion. Due to their 

'j_j ■ practical relevance, magnetohydrodynamic problems have long been the subject of intense 

Cd I cross-disciplinary research, but except for relatively simplified special cases, the rigorous 

mathematical analysis of such problems remains open. 

In magnetohydrodynamic flows, magnetic fields can induce currents in a moving con- 
ductive fiuid, which create forces on the fluid, and also change the magnetic field itself. 
There is a complex interaction between the magnetic and fiuid dynamic phenomena, and 
both hydrodynamic and electrodynamic effects have to be considered. The set of equa- 
tions which describe compressible viscous magnetohydrodynamics are a combination of the 
compressible Navier-Stokes equations of fluid dynamics and Maxwell's equations of elec- 
tromagnetism. In this paper, we consider the full system of partial differential equations 
for the three-dimensional viscous compressible magnetohydrodynamic flows in the Eulerian 
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coordinates ([HI [20]): 

Pt + div(pu) == 0, (1.1a) 

{pu)t + div (pu (g) u) + Vp = (V X H) X H + div*, (1-lb) 

£t + div{u{£' +p)) =div((ux H) x H + i/H x (V x H) + u* + kV6I), (1.1c) 

Hf- V X (ux H) = -V X (i/V X H), divH = 0, (l.ld) 

where p denotes the density, u € K'^ the velocity, H G M'' the magnetic field, and 9 the 
temperature; ^ is the viscous stress tensor given by 

* = p{Vu + Vu^) + Adivul, 

and £ is the total energy given by 

£^p(e+ i|up^ + i|Hp and £'=p(^e+ i|up 

with e the internal energy, |/o|up the kinetic energy, and ||Hp the magnetic energy. 
The equations of state p = p{p,6), e = e(p, 0) relate the pressure p and the internal 
energy e to the density and the temperature of the flow; I is the 3x3 identity matrix, 
and Vu"^ is the transpose of the matrix Vu. The viscosity coefficients A, p of the flow 
satisfy 2/i + 3A > and p > Q; v > {) is the magnetic diffusivity acting as a magnetic 
diffusion coefflcient of the magnetic field, k > is the heat conductivity. Equations (|l.lap . 
(jl.lbp . p.lcp describe the conservation of mass, momentum, and energy, respectively. It 
is well-known that the electromagnetic fields are governed by the Maxwell's equations. 
In magnetohydrodynamics, the displacement current can be neglected ([ini H^)- As a 
consequence, the equation (jl.ldp is called the induction equation, and the electric field can 
be written in terms of the magnetic field H and the velocity u, 

E = i/VxH-uxH. 

Although the electric field E does not appear in the MHD system (jl.ip , it is indeed induced 
according to the above relation by the moving conductive flow in the magnetic field. 

There have been a lot of studies on magnetohydrodynamics by physicists and math- 
ematicians because of its physical importance, complexity, rich phenomena, and mathe- 
matical challenges; see [31 HI [HI [Zl dHl (HI IHl US US] and the references cited therein. In 
particular, the one-dimensional problem has been studied in many papers, for examples, 
[21 m [T] [T31 [TS] [531 [13 and so on. However, many fundamental problems for MHD are 
still open. For example, even for the one-dimensional case, the global existence of classical 
solutions to the full perfect MHD equations with large data remains unsolved when all the 
viscosity, heat conductivity, and magnetic diffusivity coefficients are constant, although 
the corresponding problem for the Navier-Stokes equations was solved in |17| long time 
ago. The reason is that the presence of the magnetic field and its interaction with the 
hydrodynamic motion in the MHD flow of large oscillation cause serious difficulties. In 
this paper we consider the global weak solution to the three-dimensional MHD problem 
with large data, and investigate the fundamental problem of global existence. 

More precisely, we study the initial-boundary value problem of (jl.ip in a bounded spatial 
domain fJ C M'^ with the initial data: 

(p,pu,H, 6')|t=o = (po,mo,Ho,6'o)(2;), xeVl, (1.2) 

and the no-slip boundary conditions on the velocity and the magnetic field, and the ther- 
mally insulated boundary condition on the heat fiux q = —kV9: 

u|9n = 0, H|aQ = 0, q\an = Q- (1.3) 
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The aim of this paper is to construct the solution of the initial-boundary value problem of 
(|l.ip - (|1.3p and establish the global existence theory of variational weak solutions. In Hu- 
Wang [TB] , we studied global weak solutions to the initial-boundary value problem of the 
isentropic case for the three-dimensional MHD flow, while in this paper we study the full 
nonisentropic case. We are interested in the case that the viscosity and heat conductivity 
coefficients fi — fi{0), X = A(0), k = k{9) are positive functions of the temperature 6; and 
the magnetic diffusivity coefficient i^ > is assumed to be a constant in order to avoid 
unnecessary technical details. As for the pressure p = p{p, 9), it will be determined through 
a general constitutive equation: 

p = p{p,e)^Pe{p) + epg{p) (1.4) 

for certain functions pe, pg G C[0, oo) n C^(0, oo). The basic principles of classical thermo- 
dynamics imply that the internal energy e and pressure p are interrelated through Maxwell's 
relationship: 

de 1 / dp\ de dQ 

^p = 7V'^0)^ 00 = ^ = '"'^^^' 
where Cv{0) denotes the specific heat and Q = Q{9) is a function of 6*. Thus, the constitutive 
relation (|1.4p implies that the internal energy e can be decomposed as a sum: 

e{p,9)=Peip)+Q{0), (1.5) 

where 

Peip) = I' ^dt Q[9) = ^' c^{Od£.. 

If the flow is smooth, multiplying equation (jl.lbp by u and (jl.ldp by H, and summing 
them together, we obtain 

|(>P4w)+d.v(l,|uP„)+Vp.u ^^^j 

= div* • u + (V X H) X H • u + V X (u X H) • H - V X (i/V X H) • H. 

Subtracting p.6p from (jl.lcp . we obtain the internal energy equation: 

dt{pe) + div(pue) + (divu)p = v\V x Hp + * : Vu -^ div(KV6'), (1.7) 

using 

div(i/H X (V X H)) = v\V X Hp - V X {vV x H) • H, 

and 

div((u X H) X H) = (V X H) X H • u + V X (u X H) • H, (1.8) 

where ^ : Vu denotes the scalar product of two matrices (see (|4.10p '). Multiplying equation 
dnil by [pPe{p))' yields 

dt{pPe{p)) + div(pPe(p)u) +Pe(p)divU = 0, (1.9) 

and subtracting this equality from (II. 7|) . we get the following thermal energy equation: 

dt{pQ{9)) + diY{pQ{9)\i) - div(K(6')V6i) ^ v\^ x Hp + * : Vu - 9pe{p)AiY\i. (1.10) 

We note that in [6j, Ducomet and Feireisl studied, using the entropy method, the full 
compressible MHD equations with an additional Poisson's equation under the assumption 
that the viscosity coefficients depend on the temperature and the magnetic field, and the 
pressure behaves like the power law p^ with 7 = § for large density. We also remark that, 
for the mathematical analysis of incompressible MHD equations, we refer the reader to 
the work [TT] and the references cited therein; and for the related studies on the multi- 
dimensional compressible Navier-Stokes equations, we refer to [HI [SI IHl 112] and particularly 
[SI [5] for the nonisentropic case. In this paper, we consider compressible MHD flow with 
more general pressure, and use the thermal equation (|1.10[) as in [8] instead of the entropy 
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equation used in [6] , thus the methods of this paper differ significantly from those in [6] . 
There are several major difficulties in studying the global solutions of the initial-boundary 
value problem of (|l.ip - (|1.3p with large data, due to the interaction from the magnetic field, 
large oscillations and concentrations of solutions, and poor a priori estimates available for 
MHD. To deal with the possible density oscillation, we use the weak continuity property 
of the effective viscous fiux , first established by Lions |H] for the barotropic compressible 
Navier-Stokes system with constant viscosities (see also Feireisl [S] and Hoff [T^). More 
precisely, for fixed T > 0, assuming 

i{Pn,b{Pn),Pn) ^ {p,b{p),p) wcakly in L^{n x (0,r)), 
\(u„,H„) ^ (u,H) weakly in L^{[0,T];W^'\Q)), 

we will prove that, for some function &, 



{Pn - {X{dn) + 2/i(0„))divU„)fe(p„) ^ {p-{X{e) + 2fi{e))divu)b{p) 

weakly in L^{il, x (0,r)), where / denote a weak limit of a sequence {/nj^i in L^{n x 
(0, T)). To overcome the difficulty from the concentration in the temperature in order to 
pass to limit in approximation solutions, we use the renormalization of the thermal energy 
equation (|1.10p . More precisely, multiplying (jl.lOp by h{6) for some function h, we obtain, 

dt{pQh{0)) + div{pQk.ie)u) - AKh{e) 

^iy\v xii\^h{e) + h{9)'^ ■.vu-h{e)9pg{p)divu-h'{e)K{9)\ve\^, ^' ^ 

where 

n n 

Qh{e) = [ cuiOHOd^, Kn{e) = f niOH^d^. 

Jo Jo 

The idea of renormalization was used in Feireisl [3 [9] , and is similar to that in DiPcrna 
and Lions [5] . In addition, we also need to overcome the difficulty arising from the presence 
of the magnetic field and its coupling and interaction with the fluid variables. 

We organize the rest of this paper as follows. In Section 2, we introduce a variational 
formulation of the full compressible MHD equations, and also state the main existence 
result (Theorem 12.11) . In Section 3, we will formally derive a series of a priori estimates 
on the solution. In order to construct a sequence of approximation solutions, a thrcc-lcvcl 
approximation scheme from |16| for iscntropic MHD flow will be adopted in Section 4. 
Finally, in Section 5, our main result will be proved through a vanishing viscosity and 
vanishing artificial pressure limit passage using the weak convergence method. 



2. Variational Formulation and Main Result 

In this section, we give the definition of the variational solution to the initial-boundary 
value problem ()l.l|) - ()1.3p and state the main result. 

First we remark that, as shown later, the optimal estimates we can expect on the 
magnetic field H and the velocity u are in _ff ^-norms, which can not ensure the convergence 
of the terms |V x Hp and ^ : Vu in L^ of equation p.lOp . or even worse, in the sense of 
distributions. In other words, the compactness on the temperature does not seem to be 
sufficient to pass to the limit in the thermal energy equation. Thus, we will replace the 
thermal energy equality p.lOp by two inequalities in the sense of distributions to be in 
accordance with the second law of thermodynamics. More precisely, instead of (|1.10|1 . we 
only require that the following two inequalities hold: 

dt{pQ{0)) + div{pQ{e)u) - AK{d) > :^|V x Hp -h 1- : Vu - epg{p)divu, (2.1) 

in the sense of distributions, and 

E[p, u, , e, H] (i) < E[p, u, e, H] (0) for t > 0, (2.2) 
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with the total energy 

E[p,u,e,U] = j^^ (^p {p,{p) + Q{e) + i|u|^^ + i|H|2^ dx, 



and 



K{e) = f K{^)dc 

Jo 
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Now wc give the definition of variational solutions to the full MHD equations as follows. 

Definition 2.1. A vector {p,u,6,ii) is said to be a variational solution to the initial- 
boundary value problem (|l.ip - p.3p of the full compressible MHD equations on the time 
interval (0, T) for any fixed T > if the following conditions hold: 

• The density p > 0, the velocity u € L^([0,T];Wq' (ri)), and the magnetic field 
H e L^{[0,T];Wo'^{n)) n C{[0,T];Ll^^^{n)) satisfy the equations (fTTa . (fTTbll . 
and (jl.ldp in the sense of distributions, and 

n{pdtip + pu ■ Vif) dxdt = 0, 

for any (^ e C°°(f} x [0, T]) with (p{x, 0) = ip{x, T) = for x G n; 

• The temperature is a non-negative function satisfying 

/ / {pQ{e)dtLp + pQ{e)u ■ V(^ + Ki9)Aip) dxdt 
JoJn 

T I 

2 



JoJn 
for any ip e C^{n x (0, T)) with <^ > 0; 

• The energy inequality (j2.2p holds for a.e t E (0,T), with 

E[p,u,9,H]{0)^ f (poP,(pa) + poQ{eo) + l^-^^ + hno\A dx; 
Jn\ 2 po 2 y 

• The functions p, pu, and H satisfy the initial conditions in the following weak 
sense: 

ess lim / [p, p\i,'ii)[x,t)rj{x) dx — / (po, mo, Ho)?/dx, 
t-*o+ Jn Jn 

for any r]eV{n) := C^(f7). 

Now we are ready to state the main result of this paper. 

Theorem 2.1. Let fi C M"^ be a bounded domain of class C^^"^ for some r > 0. Suppose 
that the following conditions hold: the pressure p is given by the equation ()1.4p where pe, 
pg are C^ functions on [0, oo) and 

Ve(0) = 0, peiO) = 0, 

p'Ap)>ciip''-\ P^(p)>0 forallp>Q, (2.3) 

^Peip) < a2P'^ , Pe{p) < 0,3(1+ pi) forallp>0, 

with some constants 7 > |, ai > 0, 02 > 0, and 03 > 0; k = k{9) is a C^ function on 
[0,(X)) such that 

K(l + r) < k(6') <K(l + r), (2.4) 

for some constants a > 2, k > 0, and k > 0; the viscosity coefficients p and A are C^ 
functions of and globally Lipschitz on [0,oo) satisfying 

< p < p{e) < p, < \{e) < A, (2.5) 
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(2.7) 



for some positive constants /i, /i, X; v > is a constant; there exist two positive constants 
Cv, cZ such that 

~ < Cj^<c^(6l) < c;;; (2.6) 

and finally, the initial data satisfy 

' Pq e LT(f^), Po>0 on n, 

eoeL°°{n), 9o>e>o onn, 
^s^eiVr!), 

Hoei2(f7), divHo = inX>'(f7). 

Then, the initial-boundary value problem p.ip - (|1.3p of the full compressible MHD equations 
has a variational solution (p, u, 9, H) on Cl x (0, T) for any given T > 0, and 

p e L^{[o,nL''{n))nc{[o,T];L\n)), 

ueL^[0,T];W^^\n)), puGC([0,r];L;£,(n)), 

e e L^+\n X (o,r)), pQ^e) e l-([o,t];l1(^)), 

9pg e L2(j^ X (0,T)), pQ{e)u G L^f^ X (0,T)), 



(2.5 



Remark 2.1. In addition, the solution constructed in Theorem l2.1l will satisfy the continuity 
equation in the sense of renorniahzed solutions, that is, the integral identity 

/ / {b{p)dtf + b{p)u-Vip + {b{p)-b'{p)p)divuif)dxdt^O 

holds for any 

b <E C^[0,oo), |6'(z)z| < C2 2" for z larger than some positive zq. 
and any test function ip G C°°{[0, T] x H) with ip{x, 0) = ip{x, T) = for a; e Q. 

Remark 2.2. The growth restrictions imposed on K,p, A, and c-u may not be optimal, and 
7 > I is a necessary condition to ensure the convergence of nonlinear term pu (E) u in 
the sense of distributions. In particular, our result includes the case of constant viscosity 
coefficients, with the assumption that the coefficient A > 0. 

Remark 2.3. Our method also works for the case with nonzero external force / in the 
momentum equation. As it is obvious that in our analysis the presence of the external 
force does not add any additional difficulty, and usually can be dealt with by using classical 
Young's inequality under suitable assumptions on the integrability of the external force /. 



3. A Priori Estimates 

To prove Theorem 12. 11 we first need to obtain sufficient a priori estimates on the solu- 
tion. The total energy conservation (|2.2[) implies 



1, 



1, 



IHK dx 



< { PoPeipo) + PoQ{9a) + ^ — ^-t--|no 
Jn V 2 Po 2 



^IHnP I dx. 



(3.1) 



But the assumption (j2.3p implies that there is a positive constant c such that 

pPe{p) > cp^ , for any p > 0. 
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Thus, dSU implies that p'T, pQ(6'), ipjup and i|H|2 arc bounded in L°°{[0,T];L\n)). 
Hence, 

peL^{[0,nL'y{n)), pueL-([0,r];LT?T(f])). 
Next, in order to obtain estimates on the temperature, we introduce the entropy 

sip,e) = J' ^d( - Pe{p), with Pg{p) = I' ^dC 

If the flow is smooth and the temperature is strictly positive, then by direct calculation, 
using (|l.lap and (|1.10p . we obtain 

dtips) + div(p.su) + div (^) = i (j/|V X H|2 + * : Vu) - ^-^. (3.2) 



(3.3) 



Integrating (|3.2p . we get 

ps{x,t)dx— / ps{x,0)dx. 
n Jn 

On the other hand, assumptions (j2.3p imply, using Young's inequality, 

\pPeip)\ < c + pPeip) for some c> 0. (3.4) 



Moreover, we have 



since 



and 



p I ^2^d^ < pQ{e) for aU 6* > 0, p > 0, (3.5) 

" c.(0 






-dC < 0, if < e < 1, 



, rfC</ c„(OrfC = Q(e)-Q(i)<QW, ife>i- 
'1 s Ji 

Assuming that ps(-, 0) G L^(ri), then from (|3.3[) - (|3.5p , using the assumption (|2.4p and the 
estimates from p.ip . we get 

/ / |V6lt|2 + |vin6lpd.Tdt < C, 

which, combining the Sobolev's imbedding theorem, implies 

In 61 and 6lt are bounded in L^{[Q,T]-W'^-'^{n)). (3.6) 

Finally, we turn to the estimates on the velocity and the magnetic field. Indeed, inte- 
grating (|1.10p over n X (0,T), we get 

i-T I 

(* : Vu + j/|V X H|2)da;<it 
Tr ~ ~ (3.7) 



Jsi 



epg(p)diYudxdt+ I pQ{d){x,T)dx - / /9Q(6l)(a;,0)da;. 

Noticing that, using Holder inequality, one has 

\\Spe{p)\\L'^(n) < \\d\\L(^(n)\\P0{p)\\L^{n)- (3.8) 

Thus, from the assumption (|2.3p and estimate (|3.6p . we have 

0pe(p)€L2(r!x(O,T)). 
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The relation (|3.7p together with (|3.ip . (|3.8p . gives rise to the estimate 
(* : Vu + v\V X Hp) dxdi < C(po, Uq, 6Io, Hq). 



o^o 



The assumption (|2.5p . the fact ||V x Hj|^2 = j|VH||^2 when divH = 0, and Sobolev's 
imbedding theorem give that 

u,H are bounded in L^{[Q,T];W^''^{VL)). 

In summary, if ps(-,0) G L'^{^), the system ((TTTa|) . (jl.lbp . (jl.ldp . (fTTTUl) with the initial- 
boundary conditions (|1.3[) and our assumptions (|2.3p - (|2.7p yield the following estimates: 

' pPe{p),pQ{e) are bounded in L°°{[Q,T];L^{VL)); 
p is bounded in L°°([0, r];i'^(r2)), 

pu is bounded in L°°([0,T]; L^(17)); (3.9) 

In 6* and 6^ are bounded in L^{[Q,T];W^^^{n)); 
u,H are bounded in L'^{[Q,T\]Wl''^{n)). 



4. The Approximation Scheme and Approximation Solutions 

Similarly as Section 4 in [16j and Section 3 in ^ , wc introduce an approximate problem 
which consists of a system of regularized equations: 



pt + div(pu) = eAp, 

{pu)t + div (pu (g) u) + Vp(p, 6*) + ^Vp'^ 

9t((p + ,5)Q(0))+div(pQ(^)u) - A/Y (0) 

= (l-(5)(i/|VxH|2- 
Ht - V X (u X H) = -V X (i^V X H), 

with the initial-boundary conditions 

'Vp-n|ao = 
u|ao = 0, pu|t=o = "io,5 
V^lan = 0, 0|t=o = On^s, 
JA-\dn ~ 0, H|t=o = Ho. 



- eVu • Vp = (V X H) X H 

- 59°'+'^ 

* : Vu) -6'pe(p)divu, 
divH = 0, 



p|t=o = Pq,5, 



div*, 



(4.1) 



(4.2) 



where e and 5 are two positive parameters, /? > is a fixed constant, and n is the unit 
outer normal of 951. The initial data are chosen in such a way that 



Po,5- e C3(r!), Q<S<po,5<5^^] 
Po,5 -^ Po in i'^(fl), |{po,5 < Poll ^0, 
5 Jjj Pq ^ da; ^ 0, as (5 ^ 0; 



as (5 -^ 0: 



mo^s 



(4.3) 



Tno^ if Po,<5 > Po, 
^0, ifpo,5<po; 

00,8 &C^m, O<0<Oo^s<O; 
eo.s->eoin L^{n) as (5^0. 

Noticing that the terms :/| V x Hp and * : Vu are nonnegative, and = is a subsolution 
of the third equation in (j4.ip . we can conclude that, using the maximum principle, 9{t, x) > 
for all t e (0, T) and x e n. 

From Lemma 3.2 in |16| and Proposition 7.2 in [S], we see that the approximate problem 
(|4.ip - (|4.2p with fixed positive parameters e and 6 can be solved by means of a modified 
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Faedo-Galerkin method (cf. Chapter 7 in [9]). Thus, we state without proof the following 
result (cf. Proposition 3.1 in [5]): 

Proposition 4.1. Under the hypotheses of Theorem 12.11 and let (3 be large enough, then 
the approximate problem (|4.ip - (|4.2p has a solution (/?, u, 6',H) on 17 x (0,T) for any fixed 
T > satisfying the following properties: 

• p>0, ue L'^{[0,T];Wo'^{n)), H€L2([o,T];Wo'^(17)), the first equation in (|4ni 
is satisfied a.e on J7 x (0, T), the second and fourth equations in (|4.ip are satisfied 
in the sense of distributions on SI x (0, T) (denoted by I?'(ri x (0, T)), u and H are 
bounded in L'^{[0,T];WQ'^{n)), and, for some r > 1, 

Pt^Ape L^{n X (0, T)), pu e C([0, T]; i £^,(17)), 

S I I /+i dxdt < ci(e, S), e f [ |Vpp dxdt < ca, (4.4) 

JoJn JoJn 

where C2 is a constant independent of e. 

• The energy inequality 

+ (5 / / i/) (* : Vu + i/|V X up + 6l"+i) dxrii 
JoJn 



< / ,'11^0,51' , 1 



+ i|Hop + ^p1, + -T^P^s + iPo.s + S)Q{9oj)) dx 
/n \^ Po,5 2 7-1 p -1 ' ) 

+ ippb{p)divu dxdt 

JoJn 

holds for any ip G C°°([0,T]) satisfying 

^(.,0) = 1, V^(-,r)=.0, TPt<Ooiin, 
where, Pe{p) has been decomposed as 

Pe{p) = a2p'^ -Pb{p), 

withpb e Ci[0,(X)),p6 > 0; 
• The temperature 9 > satisfies that 

e L°'+\n X (0,T)), 6it G L'^{[Q,T];W^^'^{n)), 

and the thermal energy inequality holds in the following rcnormalized sense: 

{{p + 5)Qh{e)dtV + pQh{e)n ■ V^ + Kh{9)A^ ~ 6h{9)e"+^(p) dxdt 



iQJn 

< f f {{S- l)h{9){-^ : Vu + u\\/ X Hp) + h'{e)K{0)\\/e\'^) ipdxdt 
JoJn 

+ h{e)ep0{p)divutpdxdt- / {po^s + S)Qhido,s)'pix,0) dx 

Jo Jn Jn 

+ e f f ^p- V((Qft(0) - Qi9)h{e))^)dxdt, 
JoJn 

for any function h G C°°(R'*') satisfying 

h{0) > 0, h non-increasing on [0,oo), lim /i(^) = 0, 

4 — 'OO 

h"{Om > 2{h'm', for all ^ > 0, 



(4.5) 



(4.6) 
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and any test function ip G C'^{nx [0,T]) satisfying 

^>0, ip{-,T) = 0, V(^-n|ao-0. 
Remark 4.1. In fact, we have 

Comparing with the term :^|V x Hp in (|4.5[) . we need to pay more attention to the term 
\1/ : Vu, because the later involves temperature-dependent coefficients and thus can not be 
dealt with by the standard weak lower semi-continuity. Indeed, the hypothesis (j4.6p was 
imposed in |^ in order to make the function 

(e',Vu) ^ /i(6')* : Vu 

convex, and, consequently, weakly lower semi-continuous (the stress tensor 'I' in [9] depends 
on Vu only). In accordance with our new context, the following lemma is useful: 

Lemma 4.1. Let g{9) be a bounded, continuous and non-negative Junction from [0,oo) to 
M. Suppose that On and u„ are two sequences of functions defined on J7 and 

9n -^ a.e in fl, 

and 

u„ -^ u weakly in W^''^{fl). 

Then, 



(4.7) 



In particular, 



/ .g(6l)/i(6')|Vupda; < liminf / ,9(e'„)/i(6'„)|Vu„p da;. 
Jo "^°° Jn 

I h{9)'i : Vudx < liminf / /i(6l„)*(u„) : Vu„ dx. (4.8) 

Jn "^°° Jn 

Proof. First we show that ■\/g(6'„)Vu„ converges weakly to y^g{6)Wu in L^. Indeed, since 
y^g{9n)'S/Un IS Uniformly bounded in L^, it is enough to show 



v/g(6l„)Vu„(/)da-^ / ^/g(d)Vu(l)dx, for aU e C°°(17). (4.9) 

n Jn 



Since 9n ~^ d a.c in O, then \/g{9n)4> ^^ \/gW)4' ^-^ in ^ for all 4> G C°°{il). Thus by 
Lebesgue's dominated convergence theorem, we know that 

and (113 follows. 

Next, by virtue of Corollary 2.2 in [S], it is enough to observe that the function $ : 
(^7 '"^ h{6)^^ is convex and continuous on R+ x M. Computing the Hessian matrix of $, 
we get 

det{a^_^$} = 2^\h"{e)h{e) - 2{h'{e)f) > o, 

and 

trace{a^^^$} = C^h"{6) + 2h{e) > 0, 

provided 9 > and h satisfies (|4.6p . Thus the Hessian matrix is positively definite; therefore 
$ is convex and continuous. Then, (|4.7p is a direct application of Corollary 2.2 in [9]. 
Finally, from (j4.7p and the calculation: 

-v-i:^(f+£)^MWivul^ ,«o, 

follows. n 
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In the next two steps, in order to obtain the variational solution of the initial-boundary 
value problem (|l.ip - (|1.3p . we need to take the vanishing limits of the artificial viscosity 
e — » and artificial pressure coefficient 5 — > in the approximate solutions of (|4.ip - 
(|4.2p . As seen in l8l[9l[T6], the techniques used in those two procedures are rather similar. 
Moreover, in some sense, the step of taking e ^- is much easier than the step of taking 
(5^0 due to the higher integrability of p. Hence we will omit the step of taking e ^ 
(readers can refer to Section 5 in [TB] or Section 4 in [S] ) , and focus on the step of taking 
J ^ 0. Thus, we state without proof the result as e — ^ as follows. 

Proposition 4.2. Let /3 > be large enough and 5 > be fixed, then the initial-boundary 
value problem (|l.ip - (jl.3p for full compressible MHD equations admits an approximate 
solution (/5, u, 6',H) with parameter 5 (as the limit of the solutions to (|4.ip - (|4.2p when 
e — + 0) in the following sense: 

• The density p is a non-negative function, and 

pGC([0,r];L(^,,,(f7)), 

satisfying the initial condition in (|4.3p . The velocity u and the magnetic field 
H belong to ^^([0, T]; Wo'^(f7)). The equation pTTa|) and pTTd)) are satisfied in 
V'{ny. (0,r)) and 

5 [ [ p^+Uxdt<c{S). 
JoJn 

Moreover, p, u also solve equation (jl.lap in the sense of renormalized solutions; 

• The functions p, u, 6, H solve a modified momentum equation 

{pu}t + div {pu ® u) + V(]9(p, 61) + 5p^) = (V X H) X H -I- div^, (4.11) 

in 2?'(r2 X (0,T)). Furthermore, the momentum 

pueC([0,r];L£,(17)) 

satisfies the initial condition in (|4.3p : 

• The energy inequality 

[j^y^^) (^Hup + i|H|2 + pP,{p) + ^p'3 + {p + 5)Q{e)^ dxdt 

+ S [ [ ij {^ : Vu + iy\V X uf + 6"+'^) dxdt (4.12) 

< I fi^^^ + i|Hop+po,.Pe(po,5) + -^p[^, + (PM + <5)0(eM))dx 
Jn \^ Po,5 ^ P — t ' / 

holds for any ip £ C°°([0,T]) satisfying 

V'(-,0) = 1, ij{;T) = 0, Vt<0; 

• The temperature is a non-negative function, and 

9eL°'+\nx{0,T)), e""^^ eL^{[0,T];W^^^n)), cje(0,l], (4.13) 

satisfying the thermal energy inequality in the following renormalized sense: 

T/ 

{{p + 6)Qh{0)dt^ + pQh{d)u ■ Wip + Kh{e)Aip - 6h{e)e°'+^Lp) dxdt 



10 Jn 

< 

loJn 



f f {{5- l)h{e){^ : Vu + i/|V X H|2) + h'{9)K{e)\S/0\^) ipdxdt (4.14) 

JO Jo. 

nh{6)9p0{p)divuipdxdt - / {po,5 + 5)Qh{9o,5)ip{0) dx, 
. .1 Jn 
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for any admissible function h G C°°(R+) satisfying (|4.6p and any test function 
if e C'^{n X [0,r]) satisfying 



Remark 4.2. In Proposition |421 the second estimate in (j4.13p can be explained as follows. 
Taking 

in (|4.5p . we obtain 



CJ 



«;(es) 



-ive.i-'i^dxdt 



yo(i+^.)"+^ 

< - / j {pe+5)Qh{0e)i^tdxdt + 5 [ [ h{9,)e^+^i;dxdt 
JoJn JoJn 



nd^pg (pe ) I divUj I ■0 dxdt 



i-T 



\vp, ■ V ((g,.(0e) - Q{e,)h{e,))^) \ dxdt. 



JO 



Observing that 

i-T I 



and 



Jo 



vn+e^''^^ ^ i,dxdt<c I I , '^\'} ^A veA'^i^dxdt, 



Jo 



n^ePe (Pe ) I divUg I '0 (ia;(ii 
. .2 

< c||6'e||i2([o,T];Lfi(a))|lUe||i2([o,T];Wi'2(0))|iPe(Pe)||L«=([0,T]:L3(0)), 



and, by the hypothesis (|2.6[) . we have 



e/ / |Vp,-V[(Q,,(0,)-Q(0e)/l(ee))'/']|rfxdi 

Jo Jo 



<C£||V^||L-||Vpe||L2 (1 + 0^)^+1 



V6I, 



= C£||V/9e||L2 

where, the following property is used 



v{i + e,)^ 



L2 



a|Y7fl|2 



By Young's inequality. Remark 14. 11 and the energy inequality in Proposition 14. 1[ one has 

Thus, 
since 
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In the next Section, we shall take the limit of the other artificial term: the artificial 
pressure, as S ^ 0. 



5. The Limit of Vanishing Artificial Pressure 

In this section, we take the limit as (5 —> to eliminate the (5-dependent terms appearing 
in (|4.ip . while in the previous Section passing to the limit as e — > has been done. Denote 
by {psjUSjOsjIlsysyQ the sequence of approximate solutions obtained in Proposition 14.21 
In addition to the possible oscillation effects on density, the concentration effects on tem- 
perature is also a major issue of this section. To deal with these difficulties, we employ a 
variant of well-known Feireisl-Lions method [8l[9l[22] in our new context. 



5.1. Energy estimates. The main object in this subsection is to find sufficient a priori 
estimates. First, our choice of the initial data (|4.3p implies that, as 5 -^ 0, 

1 \mQ,s\' 



+ -|Hop + po,5Pe{po,5) + JZTJPq.s + (Po.s + S)Qi0o,5)] dx 



o \2 po,<5 

JnX 2 pq 2 J 

Hence, from the energy inequality (|4.12p . we can conclude that 

PS is bounded in L°°([0,T];L'' (17)), (5.1) 

VP5U5, H5 are bounded in i°°([0,r];L2(O)), (5.2) 

{ps + 5)Q{es) is bounded in L°^{[Q,T];L^{n)), (5.3) 

6jJ[p^^ + e2+')dxdt<c, (5.4) 

for some constant c, which is independent of S. 
Now, we take 

T ~ -t 1 

(z)(x,t) = 2_ h{d) = 



in (I4.14P to obtain 



^ \-^s- Vui- + i.|V X Hi|2) + ^^M^\s/es\^ ] dxdt 



JoJn \l + Os' " ' "' ^ (1 

i-T 



< 



2 / / -—^p0{ps)dWus dxdt + 2S [ [ 9^ dxdt 

-2 f (po.s + 6)Qi{eo,s) dx + f {ps + d)Qi{9s)iT)dx, 
Jn Jn 



(5.5) 



where 



Qi(0) = ^ ^de<O(0). 



Using the estimates (|5.3p and (|5.4p . we deduce from (|5.5p that 

' ^ ^ " '^ (*5 : VU5 -f v\V X H^n + ^ii^|V05|' ) dxdt 



< c I 1 + / / P0{ps)divus dxdt J , 
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for some constant c which is independent of S, and here the second term on right-hand side 
can be rewritten with the help of the rcnormahzcd continuity equation as 

npe{ps)divusdxdt = / / dtipsPeips)) dxdt. 
,. .1 JoJn 

By the hypothesis (|2.3p and the estimate (|5.ip . one has 

/ / dtipsPe{p5))dxdt < c{n) (l+ f pj dx\ < c{n) (l+ f pj dx) < c{n,T). 
Consequently, we can conclude that 



HV05peii(r!x(o,T)), 



(1 + 0^? ' 
which, combining with the hypothesis (|2.4p . gives us that 

V ln(l + ds), yej are bounded in L^{n x (0, T)). 

Thus, combining with Sobolev's imbedding theorem, we obtain that 

ln(l + 6I5), e] are bounded in L^{[Q,T]-W^^'^{n)). (5.6) 

Moreover, in view of the hypothesis p.3p . we get 

Ospeips) is bounded in L^{Vl x (0, T)). (5.7) 

With ()5.7p in hand, we can repeat the same procedure as above, taking now 

HO)^ (i + g)c. ' ^G(o,i)> ^{x,t) = '^^, 

and finally we can get 



fl ((IT^^*^- ^ ^"^- + '^'^ ^ ^^-'^^ + -(TW^I^^^'O '^'' - ^' 



(5.8) 



for some constant c which is independent of 5. 

Letting u; ^ and using the monotone convergence theorem, we deduce that 

U5, H5 are bounded in L^{[Q,T];Wq'^{VL)). (5.9) 

Moreover, from ()5.8p . wc have 

(1 + 615)''^^ is bounded in i2([o,r];VF^'2(17)), for any w e (0, 1]. (5.10) 

In particular, this implies that 

es is bounded in L^{[Q,T]-L^''+^{n)). (5.11) 

Using Holder inequality, we have 

iiC^ik^P) <ii C^iiL3(i.)ii4ii^i(^) <ii Ch\L^iD)\m\i^j,) 

for any D C il, which, together with (|5.10[) . (|5.3p and the hypothesis (|2.6p . yields 

9'^^^ dxdt<c{d), 

{ps>d} 

for any d > 0. In particular, 

/ 0'^+Uxdt < c{d) (5.12) 

for any d > 0. 
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5.2. Temperature estimates. In order to pass to the limit in the term K{9), our aim in 
this subsection is to derive uniform estimates on dg in L°'^^{n x (0,T)). To this end, we 
wiU follow the argument in [5] . 
To begin with, we have 



l{ps>d} 

where Ms denotes the total mass 



r M 

/ PS dx > Ms - d\n\ > — - d\n\, 



Ms ^ Ps dx, 
In 



independent of i e [0, T], and 

On the other hand. Holder inequality yields 



M = / padx> 0. 
n 



. 7-1 



psdxdt < \\ps\\L''{n)\{p5 > d}\ -< . 

{P5>d] 

Consequently, there exists a function K = K{d) independent of (5 > such that 

M 



\{p& > d}\ > A(d) > for aU t e [0, T], if < d < 



2|r2| 

Fix < d < ^Kt and choose a function b G C°°{R) such that 

b is non-increasing; b{z) = for 2 < d, b{z) = — 1 if z > 2d. 

For each t S [0, T], let 77 ~ rj{t) be the unique solution of the Neumann problem: 

I Ar; = b{ps{t)) -j^!^ b{ps{t)) dx in 17, ^^ ^^^ 

[V77 -11100 = 0, J^r]dx = 0, 

where 77 = ri{t) is a function of spatial variable x € fl with i as a parameter. Since the 
right-hand side of (|5.13p has a bound which is independent of (5, there is a constant r] such 
that 

Tj = r]{t) > 2 for alH > and t e [0, T]. 
Now, we take 

^ix,t) = 4^{t){ij-ri), 0<V'<1, ^P€V{0,T), 
as a test function in (|4.14p to obtain 

JJ^^ K\{es) (^b{ps{t)) - p ^^ Hpsit))] i^dxdt 



<2||77||l.o(ox(o,t)) / / {6h{es)e]+''+h{es)9spe{ps)\diYUs\)dxdt 
jQJn 

+ II'^'?IIl'°(Ox(o,t)) / / P5Qhi0s)\us\dxdt 

+ [ [ {{ps + S)Qh{Os){ri - yi)i>t - [ps + 5)Qu{es)^dtri) dxdt. 
JoJn 

Next, taking 



(5.14) 
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for < uj < 1, letting tj — > 0, using Lebesgue's dominated convergence theorem and the 
estimates (|5.ip . (|5.3p . (|5.4p . and (|5.7p . we obtain 



(5.15) 







11 

JoJn 


K{es) Upsit)) - 

V 


1 

\n\ 


/ h{ps{t))\ dxdt 






< c 
\ 


1+ / fiPs + S)9s\dt 
JoJn 


-rj\ dxdt 
J 


On 


the other 


hand, 












/ / KiOs) (b{ps{t)) ^ 


■L 


h{ps{t))\ dxdt 






= / K{es) (b{ps{t)) 

J{ps<d} \ 


1 


-Jjipsit))^ dxd 



Ki9s){b{ps{t))~— / b{psit))\ dxdt. 
{ps>d} \ l"l Jn J 

where, by virtue of (|5.12p . the second integral on the right-hand side is bounded by a 
constant independent of i5. Furthermore, 

' h{p&) dx > -— - / b{ps) dx = ~ — - > -^TTT^ > 0. 



101 ./a ' ' - MJ{p,>2d} " ' m - \^ 



Thus, we obtain 



K{es) (b{ps{t)) -^ f b{ps{t)) ] dxdt 
{ps<d} \ l^'l Jn 

>^/ K{0s)dxdt. 

I^^l J{ps<d} 



(5.16) 



Combining (|5.15p . (|5.16p together, we get 



K{0s) dxdt <cil+ I I {ps + S)es\dtv\ dxdt ] (5.17) 

{ps<d} \ JoJn J 

with c independent of S. 

Finally, since ps is a renormalized solution of the equation p. lap , we have, 

A{dtr]) = dtb{ps) ^ TTyi dtb{ps)dx 

= {b{p5) - b'{p5)p5)diYns - diY{b{ps)\is) + 77^7 / {b'{ps)ps - b{ps))divusdx. 

I"l Jn 

Hence, 

dtTj is bounded in L^{[0,T];W^'^{n)). 

Consequently, using (|5.17p . we conclude that 



/ K{9s) dxdt < c, c independent of S 

J{ps<d} 



which, together with (|5.12p and the hypothesis (|2.4p . yields 

Os is bounded in L^+"(0 x (0,T)). (5.18) 
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5.3. Refined pressure estimates. Our goal now is to improve estimates on pressure. 
We follow step by step the argument of Section 5.1 in (TB], that is, using the Bogovskii 
operator "div~^[ln(l + p5)]" as a test function for the modified momentum equation (|4.1ip . 
Similarly to Lemma 5.1 in [16j . we have the estimate 

ip{ps, 9s) + 5p^s) ln(l + Ps) dx < c, (5.19) 

n 

and hence 

p{ps, 9s) ln(l + ps) is bounded in L^{n x (0, T)). (5.20) 

Let us define the set 

jf ^{{x,t) e (0,T) X n:ps{x,t) < fc} for fc>0 and 5 € (0,1). 

In view of ()5.1|1 and the hypothesis (|2.3|) . there exists a constant s E (0,00) such that for 
all S e (0, 1) and fc> 0, 

\mT)xn-4}\<^. 

We have the following estimate: 



nSpg dxdt = / 5pg dxdt + / 5pg dxdt 

I J Jl JOx(0,T)-jf 

<TSk'^\n\+S / Xnx{a,T)-j'^Ps dxdt 



^Ox(0,T)-jf ^g21) 



Then, by the Holder inequality in Orlicz spaces (cf. [l]) and the estimate (|5.19p . we obtain 

~lx{0,T)-J^P6 



S Xnxio.T)-Jt Ps dxdt < 6\\xn,,(o.T)-j4L^ max{l, / / M{p^)dxdt} 

JoJn JoJn 

< S (n-^ (^)\ max{l, / / 2(1 + p^^)ln(l + pf ) dxdt} 

f ( k\\^^ Ft (5.22) 

<b\N-^\-\\ max{l,(41n2)r|f7| +4/3 / / p^ln(l + p^) da;di} 



< (n-^ ( - j j max{5, (41n2)<5r|f7| +AS(3 f f pfln(l + ps) dxdt}, 

where Ljv/(r2), and Ln{Q) arc two Orlicz Spaces generated by two complementary N- 
functions 

M{s) = (1 + s)ln(l + s)-s, N{s) = e" - s - 1, 

respectively. Due to (|5.19|) , we know, if 6 < 1 , 

max{,5, (41n2)(5T|rj| + 46(3 f f pf ln(l + ps) dxdt} < c, 

JoJn 

for some c > which is independent of S. Combining (|5.2ip with (|5.22p . we obtain the 
following estimate 

rT, 



j j 5p1 dxdt < TSkf^ini + c (n-^ ( - j 



where c does not depend on S and k. Consequently 



Van sup 



0^0 



6p^ dxdt 



<c(A^-(^))". (5.23) 
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The right-hand side of (|5.23p tends to zero as fc ^ oo. Thus, we have 

hm / / 5pl dxdt — 0, 
■^■^o Jo Jo. 

which yields 

(5pf -> in V'{n X (0, T)). (5.24) 

5.4. Strong convergence of the temperature. Since ps, u^ satisfy the continuity equa- 
tion (jl.lap . then 

ps^pmC{[Q,T]-Ll^^^{^)), (5.25) 

VLs^n weakly in L^{[Q,T];Wq'^{VL)), (5.26) 

and thus 

psns -^ pu weakly-* in L°°{[Q,T]-L^{n)), (5.27) 

where the limit functions p > 0, u satisfy the continuity equation (|l.lap in I>'(fi x (0, T)). 
Similarly, since ps, us, &s, H^ satisfy the momentum equation (14. lip , we have 

psus^puinC{[0,nL^^,{n)) 

and 

psus (g> us -^ pu <S) u weakly in L'^{[0,T];L^^^ (n)). 
From the hypothesis (|2.3p . equation (|l.ldp . and the estimates (|5.ip . (|5.9p . we can assume 

Peips) -^ Mp) weakly in i°°([0, T];L''{n)), 

Hs ^ H weakly in L^i[0,T];W^^'in)) n C([0,T]; LLafc(")), 
with divH = in 2?'(ri x (0, T)). Hence, p, pu, H satisfy the initial data p.2p . Due to the 
estimate (|5.10p . we can also assume 

9s ^e weakly in L^{[0,T];W^-'^{n)), 

with 6* > in V'{n x (0,T)), since 

\yOs\ <il + 9s)^"^\Ves\, forcje(0,l). 

Thus 

Ospeips) -^ OMp) weakly in L\[0,T];L^{n)), (5.28) 

(V X Hs) X Ha ^ (V x H) x H in V'{n x (0, T)). (5.29) 

Similarly, 

V X (U5 X Hs) ^ V X (u X H) in V'{n x (0, T)), (5.30) 

V X (i/V X Hi) ^ V X (i/V X H) in V'{n x (0,T)). (5.31) 

In view of (|5.6p and the hypothesis (|2.6p . we can assume 

QhiOs) -^ Q^ weakly in L^i[0,T];W^^\n)), (5.32) 

7^(6*5) ^M(^ weakly in L2([0,r];M^i'2(n)), (5.33) 

for any M e C^[0,oo) satisfying the growth restriction 

\M'{0\<c{l + ^^~'), 
and, consequently, 

ips + S)Qh{Os) ^ pQ;:M weakly in L\[0,T];LTTe(n)), (5.34) 

since L^ ^^ W~^''^{Q), if 7 > |. 

At this stage we need a variant of the celebrated Aubin-Lions Lemma (cf. Lemma 6.3 
in [9]): 
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Lemma 5.1. Let {^nj^i be a sequence of functions such that 

{^n}^=i IS bounded in L^[[Q,T];L''{n))r\L°°{[Q,T];L^{n)), withq>-, 

D 

and assume that 

dtOn>Xn mV'iVLx (0,T)), 
where 

Xn are hounded in L1([0, T]; ^^""''■(17)) 
for certain ?n > 1, r > 1. Then {0n}^i contains a subsequence such that 

Or^^e inL'{[Q,T]-W-^^^{n)). 
With this lemma in hand, wc can show the foUowing property: 
Lemma 5.2. Let h = j^, then 

{ps + 5)Qhies) ^ pQ^ m L\[0,T];W-^^^ifl)). 
Proof. Substituting h — j^ into (|4.14|) . we get 

dt {{ps + 5)Qh{0s)) > - div{psQhiOs)ns) + div ( ^^M.S70s 

V-L + fa 

na + 1 n 

iuV'iVL X (0,r)). Since 

I^q Ps (P5 ) I divu^ I < pe {ps ) | divu^ | , 
we know that, in view of (12.31) and (15.11). 



-— pe(/35)divu5 is bounded in i^([0,T]; L''(r2)), for some r > 1, 

1 + 6*5 

and, consequently, 

— ^pe(p5-)divu5- is bounded in L^{[Q,T];W-^'''{Vl)), for all fc > 1. 
Similarly, by ([51^ 

na + l 



is bounded in LWQ,T];W-^^''{9)), for all fc > 1, 



1- 
and 

diY{psQh{es)vif,) is bounded in L1([0,T]; W^'''(f7)), for some r > 1. 
Next, by (|2.4p . we have 

:J^|V05|<c(l + 05)"-'|V0,|<c0||V0j|, if05>l, 
i + t/5 



and 

"^^^ V^^l < c(l + esT-^\Ves\ < c\V0s\, if 0s < 1, 



1 

thus, by dSH) and ([SlS]) 

div (-r^^&s) is bounded in L\[0,T];W-^^''{nj), for some r > 1. 
\1 + 0s J 

Finally, since Qh{d) < Q{0), by (|5.3|) . we deduce that 

{ps + 5)Qh{es) (^ L^{[Q,T]-L\n)). 
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Hence, combining Lemma l5. II and (|5.34p together, we have 



{ps + S)Qh{0s) ^ pQhiO) in i'([0, T]-W-^'\n)). 

n 



Lemma [S^ and (|5.33p imply 



{ps + S)Qki05)Mi05) ^ pQh{d) M{e) in L^U x (0,T)), (5.35) 

where h[9) = -p-^. On the other hand, choosing M{9) = 6*, then dQh{0) satisfies (|5.33p 
since a > 2. Hence, 



U + '5)Q,.(05)^5 ^ P^g/iW weakly in L\n x (0,T)). (5.36) 

Properties (|5.35p and (|5.36|) implies 



OQhie) = Qh{0)6, a.c. on {p > 0}, 

which yields 

es-^einL^{ny,{Q,T)). (5.37) 

Indeed, wc know that Qh{0) is strictly increasing and its derivative has upper bound, 
therefore its inverse Q^^{0) exists and has lower bound l/c^. Thus, 

/ / \Qh{e) ~ Qh{es)\'' dxdt 

<ckJ f I {Qh\Qh{e))-Ql\Qh{e5))){Qh{0)-Qh{9s))dxdt 
JnJn 

^c^ I I {d- Os){Q{e) - Q{0s)) dxdt -> 0, as (5 -^ 0. 

Therefore, 

Qh{Os) -^ QhiO), in L''{n x (0,T)), as 5 ^ 0, 
and, hence, 

Qh{Os) ^ Qh{0), a.c. in r! x (0, T), as ,5 ^ 0. 
Because Q^^{0) is continuous, we deduce that 

Os^Qn\Qh{Os))^e = Q;^\Qh{0)), a.c. in r! x (0,r), as<5^0, 

which, combining Egorov's theorem. Theorem 2.10 in [5], and the weak convergence of {Og} 
to 6* in ii(r2 X (0,T)), verifies dOT]). 

Finally, (|5.37p , together with (|5.26p and Lemma 14.11 implies 

^s = fi{es){Vus + Vuf ) + X{9s)diyusl ^ * = fJ.{e){Vu + Vu^) + A(6i)divul, (5.38) 

inX>'(rJx (0,r)), and, 

/ / h{es)-^s ■■ Vvisdxdt > I I h{0)-^ : S/udxdt. (5.39) 
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5.5. Strong convergence of the density. In this subsection, wc will adopt the technique 
in [S] to show the strong convergence of the density, specifically, 

ps^ pmL\nx {0,T)). (5.40) 

First, due to (|5.20p . Proposition 2.1 in [9] and (|5.28p . we can assume that 



p{ps,9s)~^p{p,e) weakly in L^{n x (0,T)), (5.41) 

which, together with (|5.24p - (|5.3ip . implies that 

pt + div(pu) = 0, (5.42) 



dtipu) + div(pu (g)u) + \/p{p, 61) = (V X H) X H + div*, (5.43) 

Ht- V X (ux H) = -V X (i/V X H), divH == 0, (5.44) 

inV'{VL X (o,r)). 

Similarly to Section 5.3 in [16| . we use 

(^,(x,t) = ^{t)(t)[x)A.,[Tk[ps)], V e 2?(0,T), </> e V{n), i = 1,2,3, 

as test functions for the modified momentum balance equation (j4.1ip . where At can be 
expressed by their Fourier symbol as 



i^[- 



z = l,2,3. 



and Tk{p) are cut-off functions, 

Tk{p) = mm{p,k}, k > 1. 
A lengthy but straightforward computation shows that 

/ / V0 ({p{ps,0s) + Sp^g - X{es)divus)Tkips) - 2p{9s)^n,^,[Tk{ps)]] dxdt 

= 11 ^d,A [X{0s)diyus - pips, 9s) - Jpf ) A[Tk{ps)] dxdt 

+ J J ^p{es) (1^ + 1^ ) d^MATkips)] dxdt 

- J J (j)psul(i'tA[Tkip5)] + M^[iTkips) - Uips)ps)diyus]yxdt (5-45) 

nippsululdxj(pAi[Tk{ps)] dxdt 
+ / / ipug (Tk{ps)Tl.ij[psul] - (j)psulTZij[Tk{ps)]] dxdt 

- [ [ V0(V X Hs) X Us ■ Alnips)] dxdt, 
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where the operators TZtj = dxAi[v] and the summation convention is used to simpHfy 
notations. On the other hand, following the arguments of Section 5.3 in [16j . one has 



oJn 



^<P (p{p,0) - A(0)divu) Tfc(p) - 2fi{e)-^TZ^,j[Tkip)] dxdt 



^dx^cjy (A(0)divu - p{p, 0)) AATk{p)] 



dxdt 



10 Jn 

rT 



V^m(^) ( ^ + ^\dxAAATk{p)] dxdt 



loJfi 



10 Jn 
rT 



dxj dxi J 

<Ppu' (^^tA^[Tk{p)] + M'A[Tk{p) -Tl{p)p)AiYu\^ dxdt 



(5.46) 



10 Jn 

rT, 



^pu'u^dx^(j)A,[Tk{p)]dxdt 



+ i^u' {Tk{p)'R,.,[c^pu^]- (j)pu''R,,,[Tk{p)] 

JoJn ^ 

- f [ ^j(t){\I X H) X H • A[r^] dxdt. 
JoJn 



dxdt 



Now, following the argument in Section 5.3 in |16| . the Div-Curl Lemma can be used in 
order to show that the right-hand side of (|5.45p converges to that of (|5.46p . that is 



limyy V0 ({p{P6,0s) - X{es)dWus)Tk{ps) - 2p{0s)^n^,,[n{ps)] ] dxdt 



loJn 
Noting that 



dui 



(5.47) 



U ipiP,o) - x{e)diYu)Tk{p) - 2p{e)—n^^j[Tk{p)] dxdt. 



I I vp{0s)^'Jl.,,j[Tk{p5)]dxdt 



loJn 

rT 



/v, , 



-"(^S 



(pp{ds)TZi, 



du\ 



dXn 



Tk{ps)dxdt 



ipp{95)divusTk{ps) dxdt, 



oJn 



for any ip e Vl^l x (0, T)), we have, using also (|5.47|) . 



lim 

5^0 



oJn 



(p{{p{ps,es) ~ {X{es) + 2p{es))divu5)Tk{p6))dxdt 



(5.48) 



Jn 



if [p{p,6) - {X{e) + 2p{e))diyu\ Tkip) dxdt, 



since Lemma 4.2 in [8] and the strong convergence of the temperature give 



oJn 



7^^ 



(pp{0s) 



9u| 
dx-j 



(pp{0s)TZi, 



du\ 



dxi 



oJn 



n. 



sj 



9u* 



ipp{.o)n, 



■J J 



Tk{ps)dxdt 
du 



dxi 



Tk{p) dxdt. 
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As in Section 5.4 of [16], wc can conclude from (|5.48p that there exists a constant c inde- 
pendent of k such that 

Hmsup \\Tk{ps) - T'fc(p)|lL-'+i((o,T)xn) < c 

This impHes, in particular, that the limit functions p, u satisfy the continuity equation 
(|l.la[) in the sense of renormalized solutions (cf. Lemma 5.4 in [l6]V 

Finally, following the argument as in Section 5.6 in [16], (|5.40p is verified. 

5.6. Thermal energy equation. In order to complete the proof of Theorem 1 2. 1[ we have 
to show that p, u, 9 and H satisfy the thermal energy equation (jl.lOp in the sense of 
Definition [m 

In view of (|5.18p and (|5.37p . we have, as (S — > 0, 

6*5 ^ 61 in LP{n X (0, T)), for all 1 < p < 1 + a. 

Hence, by the Lebesgue's dominated convergence theorem and the hypothesis p.4|) . we 
know, as (5 — > 0, 

Khies) ^ Kh{e) '^^ L\n X {o,T)). 

By dS^m, (jOe)) and ^^AO\i . we have, as J ^ 0, 

psQh{Os)us -^ pQh{0)u, 

PsQhiOs) ^ pQhiO), 
h{0s)9spe{ps)divus -^ h{0)9pg{p)divu, 
mV'{n X (0,T)). 

Due to the strong convergence of the temperature (|5.37p and (J5.39ll - (l5.40p . we can pass 
the limit as (5 ^ in (|4.14p to obtain 

ipQh{d)dt(p + pQh{0)vL ■ Wifi + Khid)^^) dxdt 



OJfl 



I [ {h{9){^ : Vu + i/|V X li\^)if dxdt 
JoJn 

/ / h{6)9pg{p)<\iY\iLpdxdt ' 
JoJn 



(5.49) 



PQQh{Oa)^iO)dx, 



since 



oJn 



nl+a 



h{9s)dxdt 



<6 



9]+" dxdt 



h{M)6 



ei+^dxdt. 



oJn 



{0s<M} 

which tends to zero as 6 —^ 0, because the first term on the right-hand side tends to zero 
for fixed M as (5 ^ while the second term can be made arbitrarily small by taking M 
large enough in view of (|4.6p and ()5.4p . 
Next, taking 

in (j5.49p . letting w — > 0, and using Lebesgue's donnnatcd convergence theorem, we get 

r-T 



I I {pQ{e)dtip + pQ{e)u ■ Vlp + K{e)Aip) dxdt 
JoJn 



< / / iOpeip)divu - i^|V x H|^ - ^ : Vu)ipdxdt, 



(5.50) 



for any ip £ V{Q. x (0,r)) and y' > 0, since pQh{0) < pQ{0) belongs to L^{VL x (0,r)) by 
(I5TD . (12^ . and (ISTTTI) . 
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Finally, dividing (|2.ip by 1 + 6* and using equation (jl.lap . wc get 

q 



dt{pf{0))+div{pfie)u)+divl- 

> ^(HV X HP + * : Vu) - ^ - ^P.(p)divu, 
in the sense of distributions, where 

Integrating the above inequality over Q x (0, T), we deduce 



(5.51) 



< 2 sup / pf{e)dx+ [ [ -—■pg{p)\dWu\dxdt. 
o<t<TJn JaJn J^ + f 

By Holder's inequality, the estimates (|5.ip . (|5.9p . and the hypothesis (|2.3p . one has 

n- -pe(p)|divu| dxdt < / / pe(p)|divu| rfxrfi < c. 

Similarly, by Holder's inequality, the assumption (|2.6p . and the estimates (|5.ip . (|5.18p . we 
have 

/ pf{0)dx< c pOdx <c. 
Jn Jn 

Thus, ()5.5ip and the assumption (|2.4p imply that 

Ml + 9) e L^{[0,T];W^^^{n)), 6it e L'^{[0,T];W^'^{n)). 

This completes our proof of Theorem 12.11 
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